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C*") 1. Introduction It was previously suggested that inversions are present 

in atmospheres hot enough to maintain high-opacity ab- 

Observing the occultation of an exoplanet by its host star at near- sorters (such as TiO and VO) in the upper atmospher e 

infrared wavelengths allows the direct detection of thermal emis- (Hubeny et al. 2003; Burrows et al. 2007; Fortney et al. 2008). 

sion from the planet. Since the first such measurement with the More recent results (Machalek et al. 2008; Fressin et al. 2010; 

Spitze r Space Telescope dCharbonneau et al.l 2005; D eming et all iMadhusudhan & Seager! 2010) and theory (Spie gel etall 2009), 

120051) . the occultations of numerous planets have been measured however, have challenged thi s picture. An alte rnative hypothesis 

using both Spitzer at 3.6 to 24 pm, and ground-based telescopes has been put forward bv lKnutson et alJ (120101) . who note the ap- 

at shorter wavelengths. parent correlation between stellar activity and the absence of a 

thermal inversion. The planets that lack thermal inversions orbit 
Measuring the depth of a single occultation allows us to mea- c hromospherically active stars; it is proposed that the high levels 
sure the brightness temperature of the dayside of the planet and of ultravio i et flux associated with chromospheric activity destroy 
estimate how efficiently heat is redistributed from the dayside to the h i gh . a i t i tude compounds that would otherwise cause inver- 
the mghtside. Occultation measurements made at multiple wave- sions This hypothes i s is based on a small number of systems; 
lengths enable the construction of a spectral energy distnbu- observations of more planets are required to test it. Yet another 
tion (SED) of the planet's dayside atmosphere. The atmospheric possi biii ty is that TiO and VO can be naturally low in abundance 
depth probed by observations at a given wavelength depends on if the atmosph e re i s carbon-rich, C/O > 1, a nd hence preclude 
the opacity of the atmosphere in that passband, meaning that a themal invers ions (Madhusudhan etluTOlbl) . 
spectrum can inform us of the vertical structure of the atmo- 
sphere. In particular, thermal inversions (stratospheres), where WASP-24 is a late-F star which is host to a massive (1.1 
the temperature increases with altitude, have been inferred for Mj UD ) transiting planet in a short-period (2 .34 d), prograde or- 
some planetary atmospheres on the basis of occultation data. bit (IStreet et al.ll2.010t ISimpson et al.ll201ll) . The stellar flux re- 
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ABSTRACT 

Aims. We observe occultations of WASP-24b to measure brightness temperatures and to determine whether or not its atmosphere 
exhibits a thermal inversion (stratosphere). 

Methods. We observed occultations of WASP-24b at 3.6 and 4.5 yum using the Spitzer Space Telescope. It has been suggested that 
there is a correlation between stellar activity and the presence of inversions, so we analysed existing HARPS spectra in order to 
calculate log/?^ K for WASP-24 and thus determine whether or not the star is chromospherically active. We also observed a transit of 
WASP-24b in the Stromgren u and y bands, with the CAHA 2.2-m telescope. 

Results. We measure occultation depths of 0.159 ± 0.013 per cent at 3.6 pm and 0.202 ± 0.018 per cent at 4.5 pm. The corresponding 
planetary brightness temperatures are 1974 ± 71 K and 1944 ± 85 K respectively. Atmosphere models with and without a thermal 
inversion fit the data equally well; we are unable to constrain the presence of an inversion without additional occultation measure- 
ments in the near-IR. We find logR' HK = -4.98 ±0.12, indicating that WASP-24 is not a chromospherically active star. Our global 
analysis of new and previously-published data has refined the system parameters, and we find no evidence that the orbit of WASP-24b 
is non-circular. 

Conclusions. These results emphasise the importance of complementing Spitzer measurements with observations at shorter wave- 
lengths to gain a full understanding of hot Jupiter atmospheres. 

Key words, planetary systems - Planets and satellites: atmospheres - Stars: individual: WASP-24 - Planets and satellites: individual: 
WASP-24b - Infrared: planetary systems 
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ceived by WASP-24b (2.28 x 10 9 Wm" 2 ) places t he planet well 
into thepM class proposed bv lFortnev et al.l(l2008l) . so according 
to this theory a strong inversion should be present. In this paper, 
we present Spitzer occultation observations of WASP-24b at 3.6 
and 4.5 pm and we measure log R' HK for WASP-24. 

2. New observations 

2.1. Spitzer occultation photometry 

We observed two occultations of the plan et WASP-24b by it s 
host star with the Spitzer Space Telescope dWerner et al.l l2004). 
during its warm mission, on 2011 April 16 a nd 2011 April 18 . 
We used the Infrared Array Camera (IRAC; iFazio et all 12004) 
in full-array mode (256 x 256 pixels, 1.2 " pixeU 1 ). The first 
occultation was observed at 3.6 pm (channel 1) and the second 
at 4.5 pm (channel 2). Each observation consisted of a total of 
2054 10.4-s exposures, spanning a total of 7.46 hours. 

The target was placed near the centre of the array for the 
channel 2 observations, but for the channel 1 observations it 
was placed on the pixel with x and y co-ord inates {79,55). This 
position was chosen because iTodorov et al.l (1201 Ot) found low 
intrapixel sensitivity variation during observations of HAT-P-1 
centred at {79.0,55.6}. Our observations were conducted with 
the target star on a different part of the same pixel ( centred on 

1 79.0,5 5.4}), with no positional overlap with the ITodorov et al.l 
2010) observations. A strong position-dependent flux variation 
is observed in our data(see Sec. [3). 

We used the images calibrated by the standard Spitzer 
pipeline (version SI 8. 18.0) and delivered to the community as 
Basic Calibrated Data (BCD) . We performed t he photo metry us- 
ing the method presented in lAnderson et al.l (1201 lbl) . For each 
image, we converted flux from MJy sr 1 to electrons and then 
used iraf to perform aperture photometry of the target, using 
circular apertures with radii ranging from 1.5 to 6.0 pixels. The 
apertures were centred by fitting ID Gaussians with a fixed 
FWHM (1.37 and 1.62 pixels for channels 1 and 2 respectively) 
to the marginal profiles in x and y, using a non-linear least 
squares technique. We measured the sky background in an an- 
nulus centered on the target with an inner radius of 8 pixels and 
an outer radius of 12 pixels, and subtracted it from the flux mea- 
sured within the on-source apertures. 

The flux in the sky background annulus was calculated using 
iraf's centroid algorithm, which is equivalent to calculating the 
mode of the histogram of sky background values. We note that 
warm Spitzer suffers from a systematic known as 'column pull- 
down', where a column of p ixels records anomalously high flux 
readings (see iDeming et al.ll201 71 for a discussion of systematic 
effects in warm Spitzer data). One such column pulldown associ- 
ated with a hot pixel in the 4.5 pm detector is located sufficiently 
close to our target that a it appears within our sky background an- 
nulus. This does not constitute a problem for our estimation of 
the sky background, however, since only a small fraction (a; 15 
/ 250) of pixels are affected and sigma clipping is performed as 
part of the sky background estimation algorithm, so the affected 
pixels are rejected from the calculation. 

We estimated the photometric uncertainty as the quadrature 
addition of the uncertainty in the sky background (estimated as 
the standard deviation of the flux in the sky anulus) in the on- 
source aperture, the read-out noise, and the Poisson noise of the 
total background-subtracted counts within the on-source aper- 
ture. The mid-exposure times (in the HJD (UTC) time system) 
are calculated by adding half the integration time (FRAMTIME 



value in the image headers) to the MHJD_OBS values, which 
are the start times of the DCEs (Data Collective Events). These 
times were then converted to th e BJD (TDB) system using the 
routines of Eastman et al.l d2010l) . 

We adopt the aperture radii (2.6 pixels for channel 1 and 2.4 
pixels for channel 2) that give the highest signal-to-noise, and 
the smallest residual scatter in the fitted light curves (see section 
[3]for details of our fitting procedure). The variation in both the 
channel 1 and 2 occultation depths is much smaller than 1 cr 
across a wide range in aperture radius. 

In order to account for cosmic -ray hits, we rejected any 
flux measurement that was discrepant with the median of its 20 
neighbours (a window width of 4.4 min) by more than four times 
its theoretical error bar. We also performed a rejection on target 
position. For each image and for the x and y detector coordinates 
separately, we computed the difference between the fitted target 
position and the median of its 20 neighbours. For each dataset, 
we then calculated the standard deviation, cr, of these median 
differences and rejected any points discrepant by more than 4 cr. 
We reject a total of 82 (4.0 per cent) channel 1 data points (73 on 
flux, 8 on x-position, and 1 on y-position) and 23 (1.1 per cent) 
channel 2 data points (18 on flux, 1 on x-position, and 4 on y- 
position). The post-rejection data are displayed raw and binned 
in the first and second panels respectively of Figure fl] 

2.2. CAHA 2.2-m transit photometry 

We observed a transit of WASP-24b on 2010 May 05/06 with 
the CAHA 2.2-m telescope at Calar Alto, Spain. We used the 
Bonn University Simultaneous CAmera (BUSCA) to observe si- 
multaneously in the Stromgren u, b and y passbands, and the 
SDSS z passband, with exposure times of 90 to 120 s. Further 
details of our BUSCA observation strategy may be found in 
Southwor th et alj (|20 12|). Aperture pho tometry was performed 
using the pipeline of ISouthworth et aL (2009), which us es the 
idl implementation of the daophot aper algorithm (Stetson 1987) 
provided by the astrolib library. An ensemble of comparison 
stars were used and a systematic trend in the form of a lin- 
ear function of time was removed. The b and z-band data were 
both found to exhibit an unacceptably large scatter, the cause of 
which is unknown; we use only the Stromgren u and y data in 
our subsequent analysis. These data are shown overplotted with 
our best-fitting model in Figure [2] 

3. Data analysis 

Our new Spitzer data and transit photometry were incorporated 
into a global analysis alongside the following data: a total of 
29,993 photometric data points from SuperWASP-N and WASP- 
South, spanning 20 08 March 5 to 2 010 July 26 (including the 
9750 presented by Stree t et al]|2010l who included only the data 
taken prior to 2009 A pril 28); th e three complete follow-up 
transit light curves of IStreet et all (1201 Ob from the Liverpool 
(V+R band), Faulkes North and Faulkes South Telescopes (both 
PanSTARRS-z band); the ten radial velocities (RVs) measured 
by the Nordic Optical Telescope, and the eighteen CORALIE 
RVs presented by Street et al. (2010); and the 63 HARPS ra- 
dial velocities of Simpso n et al.l d201ll) . We analysed the data us- 
ing an adaptive Markov-chain Monte Carlo (MCMC) algorithm 
(Collier Came ron et al.l 120071: iPollacco et alil2008t lEnoch et al.l 
2010; see Anderson et al.ll201 lal for a description of the current 
version of our code). 

Th e dat a that we include in o ur analysis from IStreet et all 
(1201 Ol) and ISimpson et"ai1 d201ll) were originally in the HJD 
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Fig. 1. In each of the above three plots, the 3.6 fim data (green points) are shown above the 4.5 pm data (blue points), which is 
offset in relative flux for clarity. Left: Raw Spitzer data with the best-fitting trend and occultation models superimposed. Middle: 
The same data binned in phase (A</> = 0.002 « 6.74 min) with the best-fitting trend models superimposed. Right: The binned data 
after dividing by the best-fitting trend models, and with the best-fitting occultation models superimposed. The error bar on each 
binned measurement in the panels in the middle and on the right is the standard deviation of the points within the bin. 



(UTC) and BJD (UTC) time systems. We converted these 
to the BJD (TDB) system, using the routines described by 
East man et alJ d2010h . so that there is consistency between all 
our data. In order to account for the light-travel time across the 
orbit when comparing transit and occultation photometry, we 
made a first-order adjustment to the timings of the Spitzer data, 
subtracting 36 s from the time of each observation (for compar- 
ison, we measure the mid-occultation time with a precision of 
Ills). 

The MCMC proposal parameters (which are defined in Table 
0) we used are: t c , P, (R p i/R*) 2 , tu, b, K\, T e g, [Fe/H], V^cos u>, 
sfesmio, Vv sin i cos A, Vv sin i sin A, AF3 &, and AF4.5. At each 
step in the MCMC procedure, each proposal parameter is per- 
turbed from its previous value by a small, random amount. From 
the proposal parameters, model light and RV curves are gener- 
ated and;^ 2 is calculated from their comparison with the data. A 
step is accepted if x 1 (our merit function) is lower than for the 
previous step, and a step with higher x 2 is accepted with prob- 
ability exp(-A;f 2 /2). In this way, the parameter space around 
the optimum solution is thoroughly explored. The value and un- 
certainty for each parameter are taken as the median and central 
68.3 per cent confidence interval of the parameter's marginalis ed 
posterior probability distribution respectively (e.g. lFordl l2006). 

The transit light curves were modelled using the formulation 
of iMandel & A gol (2002) and limb darkening was accounted 
for using a four-coefficient, non-linear model, employing co- 



efficients appropriate to the passband from the tabulations of 
IClareti (12000, 2004). The coefficients were determined using an 
initial in terpolation in log g* and [Fe/H] (using the values re- 
ported in lStreet et al.ll2010() . and an interpolation in T,^ at each 
MCMC step. The coefficient values corresponding to the best- 
fitting value of r«, e ff are given in Table Q] 

3.1. De-trending Spitzer data 

Both the IRAC 3.6 and 4.5 pm detectors are known to have 
an inhomogeneous intra-pixel sensitivity which results in a cor- 
relation between the measured flux of a star and its position on 
the array (e.g. Knutson et al. 2008, and references therein). This, 
combined with a small, periodic (P « 40min) pointing wob- 
ble thought to be caused by thermal cycling of a heater aboard 
the spacecraft! results in a periodic modulation of the photom- 
etry (Fig. Q}. We model this so-called 'pixel-phase' effect as a 
quadratic func tion of the x and y positions of the centre of the 
PSF, after Charbonneau et al. (2008). We include in our model 
the additional cross-term of lDesert et al.l(l2009l) as well as a term 
linear in time, such that the stellar flux relative to its weighted 
mean, df — f - /, is given by 



df — «o + a x dx + a y dy + a xy dxdy + a xx dx + a yv dy + a,dt (1) 



http://ssc.spitzer.caltech.edu/warmmission/news/21oct2010memo.pdf 
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Table 1. Limb-darkening coefficients 
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Fig. 2. Transit photometry from the CAHA 2.2-m telescope. 
The upper curve (blue circles) is the transit in the Stromgren 
u passband, and the curve below it (red squares) is the transit in 
Stromgren y, offset in flux for clarity. In each case our best fit- 
ting model is overplotted (solid line) and the residuals are shown 
beneath, offset from zero and with the same symbols as the cor- 
responding light curve. 



where dx-x-x and dy — y - y are the coordinates of the PSF 
centre relative to their weighted means, dt is the time elapsed 
since the first observation, and oq, a x , a v , a xy , a xx , a yv and a, are 
the coefficients for which we fit. 

At each MCMC step, after dividing the data by the eclipse 
model we determine the trend model c oefficients using th e sin- 
gular value decomposition routine of Pres s et all ([1992). The 
best-fitting trend models are binned and superimposed on the 
binned photometry in the middle panel of FigureQ] and the trend 
model coefficients are given in Table [2] 

3.2. Photometric noise 

We scaled the formal uncertainties of each photometric dataset 
by an individual factor in order to obtain a reduced x 2 of unity. 
This ensures that each dataset is properly weighted with respect 
to the others in the MCMC analysis and that realistic uncertain- 
ties are obtained. The errors on the transit photometry were mul- 
tiplied by factors ranging from 0.95 to 2.62, while those on the 
Spitzer channel 1 and 2 data were multiplied by 1 .070 and 1 .007 
respectively. It was not necessary to add any jitter to the RV un- 
certainties in order to obtain a reduced spectroscopic-^ 2 of unity. 
We plotted the rms of the binned residuals as a function of bin 
width to assess the presence of correlated noise in the Spitzer 
light curves (Fig. |3j. A small amount of correlated noise is de- 
tected in the 4.5^/m light curve on timescales of the order of 10 
minutes. 

Remaining correlated noise present in the occultation light 
curve can result in an underestimation of the uncertainty on 
the occultation depth. To check that our occultation depth un- 
certainties are not unrealistic, we employ the residual permuta- 
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-0.353 
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Table 2. Trend model parameters and 


coefficients 






3.6 /im 




4.5 /im 







f 


193116.39 


103427.12 


X 


79.03 


127.44 


y 


55.37 


129.41 


"o 


63.9 ± 5 


-84.8 ± 2.4 


a x 


20586 ± 18 


-756.4 ± 5 


<h 


-21433 ± 53 


-3937.5 ± 6 


axy 


-9277 ± 756 


-6791 ± 329 


a xx 


-26000 ± 1302 


7784 ± 150 


a yy 


37604 ± 891 


8236 ± 282 


a t 


-633 ± 39 


59 ±21 



See Sec l3.1l for the definitions of these parameters. 



tion or 'prayer bead' method (e.g. Gill on et aLll2007l) on each of 
our Spitzer light c urves. We impleme nt this method in the same 
way as we did in iSmith et al.1 (1201 lh . i.e. we shift the residual 
to the model to the subsequent observation before adding back 
the trend and occultation models. This process is repeated, shift- 
ing each residual two observations, and so on, resulting in a se- 
ries of light curves, each of which has the time-structure of the 
correlated noise preserved. A separate MCMC analysis was per- 
formed on each of these light curves. In order to expedite the 
process of conducting a total of 4007 analyses (1977 for channel 
1, and 2030 for channel 2), we did not fit any of the transit pho- 
tometry. Instead, we fixed the values of t c , P, {R p \/R t ) 2 , tu an d b 
to the values determined by our main MCMC analysis (the val- 
ues listed in Table|3]l. We also fixed the values of ^i , T e ff, [Fe/H], 
Vv sin i cos A, and Vv sin i sin A, so that the only proposal pa- 
rameters were those influenced by the occultation photometry, 
namely the occultation depth, AFt, g or AF4.5, and -\fecos a> and 
■\Jesin u. 

We then calculated the median of the best-fitting occulta- 
tion depths from the ensemble of new fits, as well as the un- 
certainty limits which enclose 68.3 per cent of the values around 
the median. The occultation depths from this prayer bead anal- 
ysis are AF 3 . 6 = 0.159+°-°Jjj per cent and AF 4 . 5 = 0.210+J^ per 
cent. These values and uncertainties are almost identical to those 
obtained from our global MCMC analysis, indicating that any 
residual correlated noise in our light curves has a negligible im- 
pact on our measurement of the occultation depths. 



4. Results 

4.1. Brightness temperatures 

The median values of the system parameters and their associ- 
ated lcr uncertainties from our MCMC analysis are displayed in 
Table [3] Our occultation models are shown with the de-trended 
Spitzer photometry in the third panel of Figure Q] We find occul- 
tation depths of 0.159+0.01 3 percent at 3.6 pm and 0.202+0.01 8 
per cent at 4.5 pm. We calculated corresponding brightness tem- 
peratures by equating the occultation depth with the product of 
(R p \/R*) 2 and the ratio of the planet-to -star flux dens ity (calcu- 
lated using a model stellar spectrum of Kurucz 1979) integrated 
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Table 3. WASP-24 system parameters 



Parameter 



Symbol 



Value 



Unit 



Orbital period 

Epoch of mid-transit (BJD, TDB) 

Transit duration (from first to fourth contact) 

Duration of transit ingress « duration of transit egress 

Planet-to-star area ratio 

Impact parameter 

Orbital inclination 

Semi-amplitude of the stellar reflex velocity 
Centre-of-mass velocity 
Offset between CORALIE and NOT 
Offset between CORALIE and HARPS 



Orbital eccentricity 
Argument of periastron 



Phase of mid-occultation, having accounted for light travel time 

Occultation duration 

Duration of occultation ingress « duration of occultation egress 

Relative planet-to-star flux at 3.6 jim 
Relative planet-to-star flux at 4.5 /im 
Planet brightness temperature at 3.6 pm 
Planet brightness temperature at 4.5 pm 

Sky-projected stellar rotation velocity 

Sky-projected angle between stellar spin and planetary orbit axes 

Stellar mass 
Stellar radius 
Stellar density 
Stellar surface gravity 
Stellar effective temperature 
Stellar metallicity 

Planet mass 

Planet radius 

Planet density 

Planet surface gravity 

Scaled orbital major semi-axis 

Orbital major semi-axis 

Planet equilibrium temperature 7 (full redistribution) 

Planet equilibrium temperature 7 (day side redistribution) 

Planet equilibrium temperature 7 (instant re-radiation) 



p 


2.3412132 ± 0.0000018 


d 


fc 


2455149.27535 ±0.00015 


d 


t\ 4 


0.11364 ±0.00069 


d 


f l2 ~ ?34 


0.01839 ±0.00082 


d 


(V^«) 2 


0.01 102 ±0.00013 




b 


0.669 ±0.016 




i 


83.30 ±0.30 





A', 


150.6 ±2.5 


ms~' 


Trvl 


-17803 ±0.15 


ms _1 


7CORALIE-NOT 


102.85 ±0.15 


ms~' 


Tharps-coralie 


-14.21 ±0.64 


ms~' 


e 


003r 00097 
< 0.0388 (3 0-) 




OJ 


67 +24 
-14? 





ecosu 


0.00044 + «»^ 




esinw 


0.0009™ 




<rmid-occ. 

?58 

?56 ~ hi 


0.50028^ 
0.01862™ 


d 
d 


AF 3 . 6 


0.159 ±0.013 


per cent 


AF4.5 


0.202 ± 0.018 


per cent 


Tb,3.6 


1974 ±71 


K 


7b,4.5 


1944 ± 85 


K 


v sin i 


5.86 ±0.63 


kms -1 


A 


-5.8 ±4.1 





M, 


1.154 ±0.025 


M,:, 


R, 


1.354 ±0.032 


Re 


p. 


0.465 ± 0.028 


Pe 


logg. 


4.236 ± 0.017 


(cgs) 


^..eff 


6038 ± 95 


K 


[Fe/H] 


0.099 ±0.071 


(dex) 


M pl 


1.091 ±0.025 


Mi uv 


R v \ 


1.383 ±0.039 


Riup 


P P i 


0.412 ± 0.032 


Pjup 


loggp 


3. 115 ±0.022 


(cgs) 


a/R, 


5.75 ±0.11 




a 


0.03619 ±0.00027 


AU 


7p,A=0,f=l 


1781 ± 34 


K 


7p,A=0,f=2 


2118 ±41 


K 


7p,A=0,f=8/3 


2276 ± 44 


K 



1 ^p,A=o,f = / 4 r c (f J-£ where / is the redistribution factor, with / = 1 for full redistribution, / = 2 for dayside redistribution 
and / = 8/3 for instant re-radiation (Cowan & Agol 201 1). We assumed the planet albedo to be zero, A = 0. 



over the bandpass and corrected with the Spitzer transmission 
curve. We then scaled the planetary blackbody temperature so 
that the occultation depth equalled the calculated quantity. We 
find brightness temperatures of 1974 + 71 K at 3.6 /im and 
1944 ± 85 K at 4.5 pm. These uncertainties account for only 
the uncertainty in the occultation depth, which is the dominant 
source of error. 



strain the orbital eccentricity within our MCMC analysis. We do 
not detect a non-zero eccentricity with any significance, how- 
ever. Our median value is e — 0.003 1 ^00024' an d our 3c upper 
limit to the eccentricity is 0.0388. Our best-fitting phase of mid- 
occultation, m id-occ. = 0.50028^J{J225n, corresponds to an occul- 



tation occurring 57_gf s later than expected for a circular orbit. 



4.2. Orbital eccentricity 

IStreet et all (l2010h found no evidence for a non-zero orbital ec- 
centricity of WASP-24b, and adopted a circular orbital solution. 
Given the additional constraints that occultation observations 
place on e sin w, and particularly e cos u>, we chose not to con- 



4.3. System parameters 

The system parameters resulting from our analysi s should be 
viewed as superseding previous measurements by IStreet et all 
(2010) and Simpson et al. (2011), since we include in our anal- 
ysis all of the data analysed by those authors as well as addi- 
tional, new data. We made a comparison of the parameters pre- 



A. M. S. Smith et al: Thermal emission from WASP-24b 



0.0100 



0.0010 =■' 



0.0001 
0.0100 



o 



I 0.0010 r 



0.0001 




10 
Bin width / minutes 

Fig. 3. The rms of the binned residuals for the Spitzer data at 3.6 
yum (upper) and 4.5 /mi (lower) (solid lines). The white-noise 
expectation, where the rms decreases in proportion to the square 
root of the bin size, is indicated by the dashed line in each panel. 

sented in Table [3] with the corresponding values (where given) 
from lStreet et allfeoiOl) and lSimpson et al.l(l201 lh . finding them 
to be in very good agreement. The majority of our measurements 
differ from previously published values by less than 1-cr; in no 
case is the deviation greater than 2-<x. 



5. Discussion 

5.1. Brightness and equilibrium temperatures 

The fact that the 3.6 and 4.5 pm brightness temperatures are es- 
sentially identical (Table [3j Section I4.U is suggestive that no 
strong temperature inversion exists in the atmosphere of WASP- 
24b. The observed occultation depths are well fitted by a black- 
body SED for the planet (Fig. @), with a best-fitting blackbody 
temperature of 1952 ± 52 K. 

The brightness temperatures may be compared with equi- 
librium temperatures calculated under various assumptions. The 
equilibrium temperature for a zero-albedo planet is given by 



= P TVeff 



(2) 



where / parameterises the redistribution of heat; / = 1 indi- 
cates isotropic re-radiation over the whole planet (i.e. the redis- 
tribution of heat from the dayside to the nightside is fully effi- 
cient). The case where the heat is uniformly distributed on the 
dayside of the planet, but there is no redistribution of heat to 
the nightside corresponds to / = 2. A third case, correspond- 
in g to / = | (equivalen t to e - and / — | in the nota tions 
of ICowan & Agolll201 ll and lLopez-Morales & Seagerll2007l re- 
spectively), occurs when the incident radiation is immediately 
re-radiated and there is no redistribution of heat within even the 
dayside. Because the hottest region of the dayside is most visible 
close to occultation, a deeper occultation is observed than if the 
dayside had a uniform temperature. 

The zero-albedo equilibrium temperatures for the three 
aforementioned cases are presented at the bottom of Table [3] 
Our brightness temperatures fall between the equilibrium tem- 
peratures for / = 1 and / = 2, suggesting that there is signifi- 
cant redistribution of heat to the nightside. We fi nd a dimension- 
less dayside effective temperature, the T^/Tq of ICowan & Agoll 



(201 1) (where T& is the dayside effective temperature, and To is 
the equilibrium temperature of the sub-stellar point), of 0.77 ± 
0.03. In a plot of Tj /Tq vs. the maximum-expected dayside tem- 
perature (Fig. 7 of ICowan & Agolll201 lh . WASP-24b is placed 
close to TrES-4b, amongst the main cluster of points. 



5.2. Atmosphere model 

We model the dayside thermal emission from WASP-24b us- 
in g the exoplanetary atmospheric model ing and retrieval method 
of iMadhusudhan & Seagerl (120091 l2010h . The model computes 
line-by-line radiative transfer in a plane parallel atmosphere in 
local thermodynamic equilibrium. The model assumes hydro- 
static equilibrium and global energy balance. The model atmo- 
sphere includes the major sources of opacity expected in hot 
hydrogen-dominated giant planet atmospheres, namely, molec- 
ular absorption due to H2O, CO, CH4, and C02, and continuum 
opacity due to H2-H2 collision-induce d absorption (CIA). Ou r 
molecular line-lists were obtained from Freedman et al. (2008), 
Freedman (2009, private comm.), and Rothman et al. (2005), 
and we u s e the CIA opacities from Borvs ow et al.l (11997 ) and 
iBorvsowl (120021) . The chemical composition and the temper- 
ature structure of the atmosphere are input parameters to the 
model. Given a set of atmospheric observations, we explore 
the space of chemical composition and temperature structure 
to determine the regions constrained by the data fo r a desired 
level of confidence (e.g. Madhusudha n et all 1201 lal) . However, 
given that we have only two data points at present, a unique 
fit to the composition and temperature structure is not feasible. 
Consequently, we fix the molecular mixing ratios under the as- 
sumption of chemical equilibrium w ith solar elemental abun- 
dances (e.g. iBurrows & Sharplll999l) . and explore the range of 
temperature structures required to explain the data. 

We find that both our planet-star flux ratios can be explained 
by a planetary blackbody at around 1950 K. Consequently, the 
data are consistent with an isothermal atmosphere. However, 
an isothermal temperature structure over the entire atmosphere 
is unrealist ic in radiative ly efficient atmospheres at low optical 
depth (e.g. lHansenl 12008). More realistic is an atmosphere with 
temperature decreasing outward, i.e. one without a thermal in- 
version; such an atmosphere explains the data just as well (Fig. 
|4] green curve). However, a second model with identical chem- 
ical composition, but with an inversion is an equally good fit 
to the data (Fig [4] red curve). Both models predict very effi- 
cient day-night heat redistribution; the inversion model has 50 
per cent redistribution, and the non-inversion model 40 per cent. 
Our two Spitzer data points are therefore insufficient to constrain 
the presence of a thermal inversion in the atmosphere of WASP- 
24b. The two models predict different flux ratios in the near-IR, 
so a measurement of the occultation depth at one or more near- 
IR wavelengths (e.g. J, H, K-bands) would break the degeneracy 
between the models. 



5.3. Thermal inversion 



iKnutson et alJ d2010h note a correlation between the presence / 
absence of thermal inversions in the atmospheres of hot Jupiters 
and stellar activity. Planets orbiting active stars are found to lack 
stratospheres, perhaps because the high-levels of ultraviolet flux 
associated with chromospherically active stars destroys the high- 
altitude sources of opacity thought to be responsible for thermal 
inversions. 
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Fig. 4. Spectral energy distribution of WASP-24b relative to that 
of its host star. The blue circles with error bars are our best -fitting 
occultation depths. The light green line is a model-atmosphere 
spectrum, based on a model which assumes solar abundances in 
thermoc hemic al equilibrium and lacks a temperature inversion, 
and the dark red line is a model with a temperature inversion. 
The band-integrated model fluxes are indicated with circles of 
the corresponding colours. The dashed black line shows a plan- 
etary black body model with a temperature of 1900 K. Inset: 
temperature-pressure profiles for our models. 



We determined the log R' HK activity index dNoves et al.l 19841) 
of WASP-24 by measuring the weak emission in the cores of 
the Ca ii H+K lines , using the HARPS spectra obta i ned b y 
Simp son et"ai1 (1201 lh and the method of lLovis et all (1201 ll) . 
Using the spectra with the highest signal-to-noise, we find 
log R' nK - -4.98 ±0.12. This places (although not conclusively) 
WASP-24 among the inactive stars (-5.5 <logR' HK < -4.9) that 
Knutson et al. (2010) finds are host to planets with inversions. 
Furthermore, we find no evidence in the long-baseline WASP 
light curves for evidence of any rotational variability (there is 
no periodic variation above 1 mmag at a confidence level of 95 
per cent). Rotational modulation is indicative of stellar chromo- 
spheric activity, and is obser ved in e.g. WASP- 19, which has 
log R' HK in the region of -4.5 dAnderson et al .11201 Id) . 

We calcul a te the model-independent, empirical metric of 
iKnutson et alJ (120101) for classifying hot Jupiter spectra, find- 
in g C — -0.005 + 0.025 p er cent /zm (adopting the notation 
of lAnderson et al.ll201 Id) . This value suggests an inversion is 
present, but is close to the boundary between inverted and non- 
inverted atmospheres ( =a -0.05 per cent p.m -1 ), On a plot of 
\ogR' HK vs. f (Fig. 5 of IKnutson et"ai]|2010l) . WASP-24b lies in 
the upper left comer of the cluster of planets classified as having 
inversions, very close to XO-lb. 

6. Conclusions 

We have detected thermal emission from WASP-24b at 3.6 
and 4.5 /mi. Our measured occultation depths correspond to 
planetary brightness temperatures which are very similar in 
the two bandpasses; a blackbody planetary model is a good 
fit to the data. Theoretical model atmospheres which have so- 
lar abundances, both with and without a thermal inversion are 
also a good fit (Fig. @). We are unable to constrain the pres- 
ence of the therm al inversion in the a tmosphere predicted by 
both the TiO/VO dFortnev et all 120081) and the stellar activity 
dKnutson et alJJ 2010) hypotheses. 



With log R' HK = 4.98+0. 12, WASP-24 appears to be chromo- 
spherically inactive, however the relatively large uncertainty on 
this value means that the star is close to the boundary between 
active and inactive stars, although nothing else about the star 
suggests that it has an active chromosphere. This system would 
therefore present a challenge to the stellar activity hypothesis if 
near-IR occultation measurements were found to indicate that 
the atmosphere lacks an inversion. 

We find the orbit of WASP-24b to be consistent with a cir- 
cular orbit; we place a 3<x upper limit of 0.039 to the orbital 
eccentricity. 
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